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Empirical evidence is mounting that good urban design fosters the formation of social fabric. Existing evi¬ 
dence is however limited in at least two respects. First, empirical studies have focused largely on social 
interactions taking place within the residential neighborhood, while leaving social encounters near the 
workplace unconsidered. Second, while various studies have examined the impact of the built environ¬ 
ment on realized social behavior, there is as yet no empirical research on the potential for having social 
contact. A deeper understanding of the geography of social interaction potential is nonetheless important, 
for it is individuals' social opportunities rather than their preferences and actual choices that are most 
directly amenable to policy intervention. This paper seeks to address both issues in an empirical case 
study in Flanders and Brussels (Belgium). An exploratory spatial analysis is conducted to uncover spatial 
trends in the potential for social interaction in order to better understand the role of urban spatial struc¬ 
ture in the production of social interaction potential. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

It is a truism that social cohesion and community interaction lie 
at the heart of society by fostering reciprocity and trust among cit¬ 
izens (Putnam, 2000). Communities with higher levels of support¬ 
ive inter-personal interactions are likely to inspire educational 
achievement, civic engagement, economic development, respon¬ 
sive democracy, innovation and safety. At the individual level, 
being socially involved is conducive to health and well-being as 
well as the ability to find higher paying jobs. 

The need to build communities that promote social and health 
welfare has long been a central concern among planning practitio¬ 
ners and scholars (Jacobs, 1961). Empirical evidence is mounting 
that good urban design encourages the formation of social fabric. 
In particular, low-density, automobile-oriented suburban develop¬ 
ment is said to be inimical to impromptu neighborly interactions, 
while high-density, pedestrian-friendly neighborhoods with mixed 
land-uses are credited with increased levels of neighborliness and 
social vibrancy. Freeman (2001), for instance, suggests a powerful 
negative relationship between the level of car use and the number 
of social ties in a neighborhood. Leyden (2003) and Lund (2003), for 
their part, found that people living in walkable neighborhoods tend 
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to experience enhanced levels of social interaction relative to those 
living in automobile-oriented neighborhoods. 

The alleged relationships between urban design and social 
interaction are often gladly employed by New Urbanists, New 
Town advocates and other urban reformists to buttress criticism 
of urban sprawl in favor of promoting compact urban areas (San¬ 
der, 2002; Talen, 2002). Meanwhile, however, a number of scholars 
have also brought opposing evidence to the fore suggesting that 
urban sprawl in itself does not necessarily undermine social con¬ 
tact. Regressing individuals’ social interaction variables on cen¬ 
sus-tract density, among other variables, Brueckner and Largey 
(2008) observed that low-density living has a positive, rather than 
negative, effect on the propensity to engage in social activities. 
Hence, they assert that the argument of deterred social cohesion 
should not be incorporated in the panoply of critiques on urban 
sprawl. Likewise, Nguyen (2010) concluded that “compact living, 
as characterized by high population density and street accessibility 
at the county level, is unfavorable to social interaction, faith-based 
social relationships, and giving and volunteering”. 

Studies on either side of the debate, however, have at least two 
important limitations. First, they focus largely on social contact 
taking place within the residential neighborhood, while leaving 
individuals’ social encounters near the work location unconsid¬ 
ered. In other words, the implications of urban sprawl and other 
design alternatives for individuals’ commuting behavior across 
the wider metropolitan region are not accounted for. This is to be 
deemed a harmful limitation given that lengthy commutes may 
adversely affect the frequency and duration of being involved in 
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social activities (see Farber & Paez, 2009, 2011 for a detailed study 
on this matter). Furthermore, the narrow focus on neighborhood 
social ties is in sharp contrast with empirical evidence suggesting 
a tendency toward decreased intra-neighborhood and increased 
extra-neighborhood socializing (Forrest & Kearns, 2001; Guest & 
Wierzbicki, 1999). Second, the existing studies examined the im¬ 
pact of the built environment on such realized social interaction 
variables as the frequency of socializing with neighbors, member¬ 
ship in hobby-oriented clubs and attending church, but did not as¬ 
sess the potential for having social contact. Gaining insight into 
social interaction potential is nonetheless important for it is indi¬ 
viduals’ potential rather than their willingness to partake in social 
activities that is most directly amenable to policy intervention. 
This view aligns with the long-standing tradition of time geogra¬ 
phy in the social sciences which emphasizes the importance of 
constraints on (joint) activity participation (Hagerstrand, 1970; 
Lenntorp, 1978). Also, considering potential instead of observed so¬ 
cial behavior enables to circumvent the issue of self-selection (i.e. 
social people tend to prefer walkable neighborhoods, rather than 
walkable neighborhoods encouraging sociability). 

In an effort to solidify our understanding of both aspects in the 
ongoing discourse regarding the link between urban form and so¬ 
cial interaction, we have developed a method for measuring the so¬ 
cial interaction potential at the scale of metropolitan regions. The 
method is described in detail in Farber, Neutens, Miller, and Li 
(2012), but its fundamentals will be recapitulated in the next sec¬ 
tion. The method makes allowance for commuter flows and time 
budgets as well as for the land-use and transport system within 
a region. It has previously been tested with experimental data of 
synthetically constructed cities with differing land-use configura¬ 
tions and daily commute-flow characteristics. This paper presents 
an empirical application of the metrics proposed in Farber et al. 
(2012). It uses actual data from Flanders (Belgium) and Brussels 
to perform an exploratory spatial analysis of the geography of so¬ 
cial interaction potential in this region. The aim of this comprehen¬ 
sive exercise is to uncover spatial trends in the potential for social 
interaction in order to better understand the role of urban spatial 
structure in the production of social interaction potential. It should 
be emphasized, however, that this empirical study will look only at 
the space-time opportunities for having social (face-to-face) con¬ 
tact. The implications of having more such opportunities for social 
capital (Bordieu, 1986; Coleman, 1988) remain implicit and are 
merely hypothesized. The reader should thus appreciate that the 
potential for having social contact is an important, yet not the only 
prerequisite for accruing social capital. Various other structural 
and cognitive preconditions, including the quality of social re¬ 
sources (see Rostila, 2010), need to be considered to fully grasp 
the notion of social capital. 

The remainder of this paper is organized in four sections. The 
next section describes the fundamentals of the method introduced 
earlier in Farber et al. (2012). Section 3 presents the study area, 
data, assumptions and computational aspects. Results are pre¬ 
sented and discussed in Section 4. Finally, the paper concludes 
with a brief outline of the major findings as well as the avenues 
for future work. 


2. Method 

2.1. Conceptual approach 

Central to our approach (see also Farber et al., 2012) is the deri¬ 
vation of the potential for social interaction from the intersection of 
space-time prisms (Fig. 1). A space-time prism is a key concept of 
time geography and gathers all space-time points that are physi¬ 
cally accessible to an individual given one or more space-time 


anchors, the maximum velocity of physical movement and the min¬ 
imum time required for some activity (Hagerstrand, 1970). Space- 
time anchors represent key locations in an individual’s life, such 
as home and work, where activities with a high degree of space- 
time fixity are undertaken - that is, activities that are relatively 
difficult to re-schedule or re-locate in the short run. Space-time an¬ 
chors condition physical accessibility by dictating where and when 
discretionary travel and activities must start and end (Cullen & God¬ 
son, 1975; Kwan & Hong, 1998; Miller, 2005a). 

In recent years, several authors have relied on prism intersec¬ 
tions as a way to model the potential for social interaction. Recent 
accomplishments in this field include the implementation of tool¬ 
kits for querying and representing the opportunities for joint activ¬ 
ity participation and intra-household interaction (Kang & Scott, 
2008; Kwan & Lee, 2011; Neutens, Versichele, 8; Schwanen, 2010; 
Yu & Shaw, 2008), the theoretical formulation of necessary condi¬ 
tions for physical and virtual interaction (Miller, 2005b) and the 
development of an analytical framework for measuring joint acces¬ 
sibility (Fang, Tu, Li, & Li, 2011; Neutens, Schwanen, & Miller, 2010; 
Neutens, Schwanen, Witlox, & De Maeyer, 2008). 

In line with these studies (Fang et al., 2011; Neutens, Schwanen, 
& Miller, 2010; Neutens et al., 2008), we will express the potential 
for social interaction in terms of the amount of time available for 
interaction with others within a particular time budget. As com¬ 
muting constitutes the primary nexus of daily spatial mobility (Sal- 
ze et al., 2011), the subsequent theoretical and empirical 
development will concentrate on the interaction possibilities with¬ 
in a single time budget between work and home. In other words, 
individual prisms will be anchored at the home and work location 
as depicted in Fig. 1. The model can be extended to other parts of 
the day and multiple time budgets in future research, with the 
appropriate modifications. As journey-to-work travel data for an 
entire metropolitan region is typically not available at the individ¬ 
ual level, we will adopt a zonal approach and consider commuting 
flows from one zone to another. Hence, we will translate the 
microscopic rendezvous scenario between two individuals, as de¬ 
picted in Fig. 1, to a general situation where an individual with a 
particular commuting pattern may socialize with all other com¬ 
muters in a given metropolitan region. This means that we are con¬ 
cerned with the potential intersection of zone-to-zone commuting 
flows rather than prism intersections of separate individuals living 
and working at discrete anchor locations. In the next section we 
will explain more formally how these intersections can be 
computed. 
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2.2. Theoretical framework 


The amount of time available for an individual living in zone i 
and working in zone j to perform an after-work out-of-home activ¬ 
ity in zone k is expressed as: 

jk _ J % ~ (tjk + hri) if tjk + t ki < t bjj ^ 

11 | 0 otherwise 


where t bj . = t h .. - t w . is the time available for a discretionary activity 
after leaving work and before arriving home; t jk is the time it takes 
to travel from work in zone j to zone k\ t ki is the time it takes to tra¬ 
vel from zone k back to home in zone i; t Wjj is the time the individual 
can leave work; and t hi is the time the individual must arrive back 
at home. Note that Eq. (1) ensures that time availability is con¬ 
strained to non-negative numbers. T^k can be considered the max¬ 
imum amount of time an individual can be in zone k given the 
capability constraint of physical travel through space and the 
work/home anchor constraints. 

Let T k ? = t Wij + tj k be the start time of the time window of T| and 
let T k t = t hjj - t ki be the last moment of time an individual may be 
in zone k before needing to return home. Likewise, T q s r = t Wqr + t rk 
and T q e r = t hqr - t kq denote the start and end times of the time win¬ 
dow in zone k for an individual living in zone q and working in 
zone r. 

It then follows that the amount of time these two individuals 
can interact in zone k is: 


= f min(T^T^) -max(T^,T^ r ) if min(T^,T^ r ) >max(7f ,7^) 

1 0 otherwise 


( 2 ) 


Again, Eq. (2) ensures that the time duration for a joint activity is 
non-negative. Furthermore, given a finite set K of discrete potential 
interaction zones, the total potential time available for after-work 
interaction between these two individuals can simply be expressed 
as: 


E 4 r 


( 3 ) 


Now suppose there is a metropolitan region with zones and a work¬ 
ing population R. Let |R| = Ei/s where r,j is the number of workers 
living in zone i and working in zone j. Additionally, let P q = p be the 
share of workers in the metropolitan region that travel from zone i 
to zone j. 

Then, 

= ££££(wvv) ( 4 ) 

i j q r 


constitutes a metric of social interaction potential for metropolitan 
region c dependent on the spatial configuration of land-uses and 
transport infrastructure, and the realized commuting flows of the 
region. M c can be interpreted as the aggregate of potential social 
interaction time between all pairs of commuters in the city, 
weighted by the likelihood of occurrence of each type of commute. 

The metric M c attains a maximum when travel velocities are 
infinite, or when all residential, employment, and social activity 
locations are exactly collocated and no travel is needed in the sys¬ 
tem. This maximum is given by 

= Kc££££PiAr [min(t hli , t hqr ) - max(t w .., t w , r )l (5) 

i j q r 

where K c is the number of zones defining metropolitan region c . 

The metric attains its minimum of 0 when the required time to 
reach all social activity locations is greater than the allowed after¬ 
work time budgets (i.e. Vi,j,k : tj k + t ki > t b ). 


In order to make absolute amounts of social interaction 
opportunity of different regions commensurable, normalization 
by the number of social interaction zones identified in each 
region is required. Dividing M c by K c modifies the interpretation 
of the metric to be the expected amount of social interaction time 
at any location between any two commuters in the considered 
urban region.M c can be decomposed into four meaningful met¬ 
rics which will be employed to examine the social potential of 
zones internal to a given metropolitan region These metrics 
include: 


= £££(W.) (6) 

j q r 

M a = eeeowjw) ( 7 ) 

i Q r 

Mdj = ££(« p ,r) (8) 

q r 

^=££££( 4 *) o) 

i j q r 


Notice that M c = E.Mci = Ej M cj = EiEj M dj = Ek M ck■ The first 
metric will be used to determine the social potential of living in 
zone i of urban region c. The second for determining the social po¬ 
tential of having a job in zone j. The third for determining the social 
potential associated with a given daily commute from zone i to j. 
And the fourth will be used to determine the amount of potential 
social interaction demand in zone k. 

For a more elaborate discussion of the scope and interpretation 
of the proposed metrics as well as their broader applications, the 
reader is referred to Farber et al. (2012). 


2.3. Example 


As an illustration, we will show how M cq (the metric with the 
least iterations, see Eq. (8)) can be calculated in practice using a 
simple example of a fictitious metropolitan area containing 100 
workers and consisting of three zones A, B and C. Commuting flows 
as well as travel times within and between the zones are presented 
in Fig. 2. For example, zone B attracts 20 workers from zone A and 
25 from zone C, who need to travel 12 and 20 min to their work 
place, respectively. Zone B also employs local workers from within 
the zone who on average travel 7 min to their workplace. It is 
assumed that the time budgets of all individuals within the 
metropolitan area equal 90 min, starting at 5 PM and ending at 
6.30 PM. 

The social interaction potential associated with a daily commut¬ 
ing flow from zone A to B can then be calculated as follows: 


McAB = y^y^ (7/tflqrP/'BfV) 
q r 

= TabmPabPm + TababPabPab + TabacPabPac + TabbaPabPba 
+ TabbbPabPbb + TabbcPabPbc + TabcaPabPca 
+ TabcbPabPcb + TabccPabPcc 

_ T 20 10 20 20 20 5 

ABAA 100 100^ ^100 100 + ^iooioo 

20 5 20 15 20 5 

+ M“ 10 0 100 + ABBB 100 100 + ABBC 100 100 

20 10 20 25 20 5 

+ A fl « 100 ioo+ * M 100 100 + ABCC 100100 


( 10 ) 


In this equation Tabaa should be calculated according to Eq. (3) as 
follows: 
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Fig. 2. Calculation example. 


r _ \ fk _ -Wl | 'rli , J'L 

l ABM — / ,>ABM — 1 ABM ' ‘ ABM ' 1 ABM 
keK 

= (min(6.24 PM,6.24 PM) - max(5.13 PM.5.06 PM)) 

+ (min(6.17 PM,6.17 PM) 

- max(5.07 PM.5.07 PM)) 

+ (min(6.15 PM,6.15 PM) 

- max(5.22 PM,5.17 PM)) 

= 71 min + 65 min + 53 min = 189 min (11) 

Eq. (11) denotes the total amount of time that a commuter from 
zone A to B would be able to interact in the entire metropolitan area 
with a commuter who works and lives in zone A. These two com¬ 
muters can spend together 71 min, 65 min and 53 min in zone A, 
B and C, respectively. The total potential time available for after¬ 
work interaction between other commuter type pairs ( T ABAB , T A bac, 
Tabb At etc.) can be calculated analogously. So doing, it then follows 
that the social interaction potential associated with a daily com¬ 
mute from zone A to B is: 

M cA b = 189 min • 0.02 + 194 min ■ 0.04 + 179 min ■ 0.01 

+ 178 min ■ 0.01 + 183 min ■ 0.0375 + 168 min ■ 0.01 
+ 169 min ■ 0.02 + 174 min ■ 0.05 + 159 min ■ 0.01 
= 37,32 min (12) 

Eq. (12) denotes the contribution associated with the commuting 
flow from A to B to the overall social interaction potential in metro¬ 
politan region c. The other metrics of social interaction (i.e. M ci , M CJ 
and M ck ) can be calculated in a similar fashion. Importantly, one can 
already see that even for this simple example consisting of three 
zones, the computation gives rise to a large number of iterations. 
We will explain how we have dealt with this computational burden 
in the next section, after discussing the characteristics of the study 
area and the data. 


3. Data 

3.1. Study area and context 

The study area consists of Flanders and Brussels Capital Region 
(hereafter simply referred to as Brussels). Flanders is the Dutch¬ 
speaking, northern half of Belgium (Fig. 3). The area is far more 
densely populated than its southern, French-speaking counterpart 
(Wallonia) (Fig. 4). Flanders is approximately 13,500 km 2 in size 
and has over 6 million inhabitants (i.e. about 60% of the Belgian 
population). Flanders is one of the six Central and Capital City Re¬ 
gions (CCC regions) of international economic importance in 
North-West Europe, alongside London, Paris, the Frankfurt Rhine- 
Main region, the Rhine-Ruhr region and Randstad Holland. Being 
one of the most densely urbanized areas in the world, Flanders 
forms a challenging study area to implement the approach out¬ 
lined in Section 2. While Brussels (160 km 2 ) is part of neither Flan¬ 
ders nor Wallonia, it will be included in the analysis because of its 
huge importance in the context of commuting flows. 

According to the Spatial Structure Plan of Flanders (RSV, 2004), 
the study area includes three metropolitan regions: Brussels, An¬ 
twerp and Ghent. The major nexus of economic development is 
formed in and along the connections between these metropolitan 
regions and is known as the Flemish Diamond. Brussels (over 1 mil¬ 
lion inhabitants) is the largest employment center in Belgium and is 
mainly characterized by service industries. The service industry is 
also strongly represented in the center of Antwerp (approximately 
500,000 inhabitants). Furthermore, Antwerp accommodates one 
of the largest ports in Europe and constitutes an important petro¬ 
chemical node, with seven of the ten largest chemical companies 
in the world having a production site in this area. Ghent (approxi¬ 
mately 250,000 inhabitants) is a medium-sized city and particularly 
known as a research and development center. Beyond the Flemish 
Diamond, there are a number of other significant smaller urban re¬ 
gions which, together with Antwerp and Ghent, are termed ‘central 
cities’ (see Fig. 4), referring to the fact that they serve as a central 
node for their surrounding areas in terms of employment, educa¬ 
tion, services, culture and recreation. The two most important 
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Fig. 4. Population density in Belgium. 


non-metropolitan urban regions in Flanders are Bruges (120,000 
inhabitants) and Leuven (95,000 inhabitants). The former is well- 
known as a tourist destination and for its port activities, in particu¬ 
lar regarding roll-on-roll-off shipping and the transferral of new 
vehicles. The latter contains the largest Dutch speaking university 
and attracts the largest share employees in the creative and knowl¬ 
edge-based industry of all municipalities in the country. In order of 


decreasing population, the other central cities include Mechelen 
(81,000), Aalst (80,000), Kortrijk (75,000), Hasselt (73,000), Sint- 
Niklaas (71,000), Oostende (69,000), Genk (65,000), Roeselare 
(57,000) and Turnhout (41,000) (Fig. 4). 

The commute areas of the various urban regions in the polycen¬ 
tric patchwork of Flanders and Brussels overlap strongly. The aver¬ 
age commuting distance in the study area amounts to 16 km (one 
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way) (Verhetsel, Thomas, Van Hecke, & Beelen, 2007). In recent 
years, there is a large and systematic rise in the average commut¬ 
ing distance: between the successive decennial censuses of 1991 
and 2001, the average commuting distance has increased by more 
than 10% (Thomas, Verhetsel, & Beelen, 2008). Importantly, almost 
one fourth of all commuters in Belgium (i.e. persons working out¬ 
side the municipality where one lives) have their destination in 
Brussels. Hence, Brussels is expected to be a focal point of after¬ 
work social interaction potential. However, while the high working 
density implies more people to socialize with, the longer 
commutes to Brussels will act counteractively by implying less 
time for social interaction. Apart from Brussels, other important 
employment centers with relatively more incoming than outgoing 
workers include Antwerp, Ghent, Liege, Zaventem (airport), Charle¬ 
roi and Leuven. Residential areas with a high negative balance be¬ 
tween the incoming and outgoing flow of workers are primarily 
located in the surroundings of the major employment centers. 
The lowest average commute times in the country can be found 
in the western part of Flanders. 

The contemporary commuting structure in Belgium has been 
shaped historically by an array of intended anti-urban policies 
and socio-economic processes that have taken place since the 
Industrial Revolution. In brief, the development of a fine-meshed 
network of railway connections and the abolishment of charges 
to enter inner cities in the second half of the 19th century drove 
many workers to the fringes of the cities that until then had had 
sharp boundaries - often materialized through medieval city walls. 
This wave of suburbanization was subsequently invigorated and 
institutionalized by a social tariff scheme for commuting by train 
and other policies that were aimed at keeping the working class 
away from the bad influences of the city and under the influence 
of the local clergyman (see De Decker, 2011 for an in-depth discus¬ 
sion of consecutive uban policies in Flanders). In the postwar era, 
increased affluence and a concomitant increase in mobility, 
changes in residential preferences and the rapid growth of the ter¬ 
tiary sector gave rise to urban flight and peri-urbanization in the 
form of detached, single-family housing on the countryside; pro¬ 
cesses which continue until today (Boussauw, Derudder, & Witlox, 
2011 ). 

3.2. Joumey-to-work flow matrix 

The journey-to-work flow matrix for Flanders is derived from 
the MultiModal Model (MMM). The MMM is a macroscopic traffic 
simulation model commissioned by the Flemish government since 
1998 (Boussauw, Neutens, & Witlox, 2011; Verhetsel et al., 2007). 
The model is calibrated on the basis of traffic count data, the 
decennial census and a region-wide activity-travel diary data set 
(OVG, ‘Onderzoek VerplaatsingsGedrag’). It simulates the number 
of trips, disaggregated according to their purpose, that take place 
during a particular time window between any pair of traffic anal¬ 
ysis zones (TAZs). For the current study, the simulated number of 
work trips made between 4 and 11 am on an average weekday 
morning has been used as a proxy for the number of workers that 
commute from one TAZ to another. 

The MMM is essentially composed of five sub-models, one for 
each province of Flanders (i.e. Flemish Brabant, Antwerp, East Flan¬ 
ders, Limburg and West Flanders). Brussels is included in the sub¬ 
model of Flemish Brabant. Every sub-model is represented by a 
CIS layer that divides the considered province into small TAZs. 
While in densely populated areas TAZs largely correspond to census 
tracts, census tracts have been aggregated into a single TAZ in more 
peripheral parts of the province. The areas surrounding the consid¬ 
ered province (i.e. Wallonia, France, The Netherlands, Germany, 
Luxembourg and the other Flemish provinces) are also represented 
in the CIS layer and divided in TAZs of which the size increases with 


increasing distance from the province under scrutiny. The CIS layers 
for every sub-model can be linked to the journey-to-work flow ma¬ 
trix by a unique identifier. The number of TAZs amounts to 1404 in 
West Flanders, 2007 in East Flanders, 2059 in Limburg, 2689 in 
Flemish Brabant and 3320 in Antwerp. To obtain a detailed map of 
the potential for social interaction for Flanders, we have applied 
the proposed method to each of the five provincial journey-to-work 
subsets separately, normalized these by the number of TAZs, and 
spatially merged the results afterwards. 

3.3. Travel time matrix 

Free-flow travel times by car along the shortest path between the 
centroids of all pairs of TAZs in each provincial sub-model have been 
estimated using ESRI®’s ArcGIS™ Network Analyst. To this end, we 
have used the TeleAtlas® MultiNet™ data set (version 2007.10). This 
data set contains a detailed topological representation of the Belgian 
road network built up by links and nodes to which various attributes 
such as turn restrictions, speed limit, and functional road classifica¬ 
tion are attached. In total, 28,491,867 shortest paths had to be cal¬ 
culated in order to construct five matrices of travel times between 
all possible pairs (i,j) of TAZs of the different partitioning systems 
used in each provincial sub-model. The provincial travel time matri¬ 
ces are quasi-symmetrical meaning that in most cases travel time 
from i to j roughly equals travel time fromj to i. 

As the above shortest path calculations generated intrazonal 
travel times equal to zero on the diagonal line of the travel time 
matrix, we have adopted the approach suggested by O’Kelly and 
Lee (2005) to approximate average intrazonal commute. The ap¬ 
proach assumes that each zone has a circular shape such that the 
average intrazonal travel distance can be estimated backwards 
from the known zonal area. When assuming an equal distribution 
of trips across zones, this intrazonal travel distance equals the ra¬ 
dius of a zone (Frost, Linneker, 8i Spence, 1998; Horner, 2002): 


where p„ is the approximate average intrazonal distance, and S,- is 
the area of zone i. Finally, pa is converted into travel time by assum¬ 
ing an average travel speed of 60 km/hour. 

3.4. Assumptions 

It is important to clearly emphasize the limitations and assump¬ 
tions underlying this study. First, as explained in Section 3.3, travel 
costs will be expressed by car travel times. While the car is cer¬ 
tainly the dominant commuting mode in Flanders and accounts 
for over 71% of all commuting trips (65% as driver and 5% as pas¬ 
senger), the share of other transport modes is also substantial: 
bicycle (13%), train (7%), bus tram and metro (4%), pedestrians 
(3%) and motorbikes (2%). Car use is relatively less dominant in 
major cities and along the coastal area where more public transit 
options are available and commuting distances are smaller 
(Verhetsel et al., 2007). Second, congestion effects and parking 
problems have not been incorporated in this study. This may lead 
to an underestimation of commuting times (and thus an overesti¬ 
mation of social interaction potential), especially for those working 
in the metropolitan areas of Brussels, Antwerp and Ghent. Third, 
given the aggregate nature of the origin-destination matrices, we 
were unable to differentiate individuals (or zones) on the basis of 
their departure/arrival time from/at work/home. Hence, a uniform 
after-work time budget of 90 min has been assumed for all com¬ 
muters in the study area. In other words, geographic proximity 
of commuting flows will be the main enabler of social interaction 
potential, while the incorporation of timing constraints is left for 


324 


T. Neutens et al./Computers, Environment and Urban Systems 41 (2013) 318-331 


future work. Fourth, our approach does not account for the fact 
that some form of social interaction can also happen during the trip 
itself (e.g. in the car, on the train), which, for some people, may 
lead to an underestimation of their social interaction potential. Fi¬ 
nally, only the after-work social interaction potential of employees 
is looked at in this paper; people who are not employed outside of 
the house have not been considered. 

The principal reason why the above assumptions have been 
(had to be) made is the lack of appropriate data to adequately mod¬ 
el these effects in Flanders. For example, estimations of working 
hours, congestion, parking times and the number of carpoolers 
are difficult to obtain in a sufficiently detailed form for the entire 
study area. Furthermore, the origin-destination matrices available 
to us can only distinguish between trip purposes not between 
transport modes. Also, the calculation of travel times by public 
transport and travel times of multi-modal trips would be very 
time-consuming and difficult given the large number of calcula¬ 
tions (almost thirty million) that have to be performed. While 
the above limitations should definitely be addressed in the future 
when finer-grained data would be available for Flanders, we be¬ 
lieve the assumptions are nonetheless reasonable in the frame¬ 
work of this first exploratory study, given that its purpose is to 
identify large-scale spatial trends rather than individual level dif¬ 
ferences in the potential for social interaction. Nevertheless, the 
reader should be cautioned that the results obtained may not nec¬ 
essarily hold under different assumptions. 

3.5. Computation 

The social interaction potential metric is a weighted average of 
all possible after-work space-time prism intersections. A region of 
N zones has N 2 possible origin/destination pairs and the same 
number of possible space-time prisms. This gives rise to N 4 poten¬ 
tial prism intersections. Since the volumes of the prisms are 
discretized over the N zones of potential social activity, the number 
of computations required is in the order of 0(N 5 ). Thus a simple re¬ 
gion with 100 zones may require 10 billion zone-based calcula¬ 
tions; increase the complexity of the zoning system to 2000 
zones, and the social interaction potential metric requires 
3.2 x 10 16 , or 32 million billion computations. In practice the num¬ 
ber of calculations is reduced by pre-screening the origin/destina¬ 


tion matrix for non-zero flows. Even with this reduction, practical 
computation of the metric necessitates a high performance com¬ 
puting infrastructure. 

The metrics of social interaction potential were computed using 
the clustered computing facilities of The Center for High Perfor¬ 
mance Computing at The University of Utah (http://chpc.utah.edu). 
The Ember Cluster used in this analysis consists of 262, 2.8 GHz, 
12-core nodes, giving rise to 3144 total cores with 6288 GB of total 
memory. This system enables fast computation of highly parallel- 
izable tasks. For example, using 768 simultaneous cores, the sys¬ 
tem computed the full cadre of metrics for Antwerp in 2 h 
35 min, for Flemish Brabant in 2 h and 30 min, for East Flanders 
in 36 min and West Flanders in 8 min. Similarly, using 960 cores, 
the metrics for Limburg were computed in 13 min. The computa¬ 
tional overhead required to partition the problem into parallel 
components and assemble the results into a unified output matrix 
is negligible. This means that the time required to compute the 
metrics is perfectly proportional to the number of processors used. 
So, for example, doubling the resources used for Antwerp would 
likely result in a total computation time of only 1 h and 20 min, 
while reducing the number of processors to one - in order to sim¬ 
ulate the use of a single-core desktop computer - would result in a 
total computation time of nearly 83 days. 

4. Results 

4.1. Social interaction potential in Flanders 

First we explore the general spatial structure of social interac¬ 
tion potential of commuters in Flanders using the four different 
metrics specified in Section 2.2 (i.e. M d , M CJ -, M ck and M C y). Given 
the extreme values obtained at the center of Brussels, the values 
of the different metrics have been categorized and visualized in 
quantiles to keep sufficient spatial variation in the figures. Fig. 5 
depicts the spatial variation in M ri (represented per zone of resi¬ 
dence i) - that is it shows how socially sustainable it is to live at 
a particular location in Flanders given the realized commuting 
flows as well as the land-use and transport system in the area. 
Fig. 5 shows a clearly articulated spatial structure for M ci . The high¬ 
est values are encountered in and around the Flemish Diamond 
formed by the mutually well-connected metropolitan areas of 
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Fig. 5. Social interaction potential of living in Flanders (M d ). 
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Fig. 6. Difference in social interaction potential in Flanders between living and having a job (M d — M cj ). 


Brussels, Antwerp and Ghent. This is not surprising given the high 
population density and the large number of people who commute 
to these cities on a daily basis. Social interaction potential is also 
relatively high for those living in the core of the central cities as 
well as in some major cities along the Belgian coast, despite their 
eccentric location in the study area. Residential areas that exhibit 
poor levels of social interaction levels, on the other hand, can be 
found in the western and north-eastern corner of Flanders, close 
to the French and Dutch border, respectively. Knowing that the 
average commuting distance in the south-western part of Flanders 
is lower than elsewhere in the study area, it appears that the influ¬ 
ence of commuting distance on social interaction potential is over¬ 
compensated by the relatively low density of both jobs and 
residents, as well as the absence of important highways in this 
area. 

To gain insight into the difference between the social interac¬ 
tion potential associated with either living or having a job in Flan¬ 
ders, Fig. 6 compares the spatial pattern of M ci to that of M CJ . One 
can observe that in Brussels and the central cities, the average so¬ 
cial interaction potential of workers is generally higher than that of 
residents - especially in the city center, industrial zone and/or port 
of these cities, where a large number of jobs are concentrated. The 
same is true for recreation areas (i.e. parks, forests, etc.) where 
fewer (if any) people live than work. 

The spatial variation of M ck (i.e. the expected potential demand 
for after-work social interaction) is represented in Fig. 7. In other 
words, Fig. 7 shows how probable it is that after-work social con¬ 
tact can take place between any two commuters at a particular 
location in Flanders. The spatial pattern resembles a more pro¬ 
nounced version of that of M ci . There is a sharp contrast in potential 
social demand between the values in Brussels and those encoun¬ 
tered in smaller cities. The high density of prism intersections in 
Brussels resonates, among others, the location strategy of various 
lobby organizations for which informal after-work social contacts 
with decision makers are paramount. Taken together, the demand 
for after-work social interaction can attain the highest levels at 
locations in the central part of the study area where the potential 
convergence of commuting flows is large. 

To further explore why some regions may have higher (lower) 
social interaction potential than can otherwise be explained by 
the density of the working population, we have regressed social 
interaction potential on a selection of potential explanatory 
variables related to commuting and urban form. These variables in¬ 


clude density of the working population 1 (dens), job-housing bal¬ 
ance 2 ( jhb ), spatial-functional diversity 3 (div), average commuting 
distance ( avgd ), (Euclidean) distance to highway ( dhigh ) and (Euclid¬ 
ean) distance to major city center 4 ( dcen ). We checked whether any 
of the variables were strongly correlated with one another (r> 0.4); 
this turned out not to be the case. However, Moran’s I and Lagrange 
Multiplier (LM) test statistics revealed that there is significant spatial 
dependence in the variables (‘spatial lag’) and/or errors terms (‘spa¬ 
tial error’) of a traditional OLS regression model. This may, among 
others, be due to the fact that the boundaries for which the data 
has been gathered do not reflect the nature of the underlying pro¬ 
cesses. Both LM-lag and LM-error test statistics are highly signifi¬ 
cant, with the latter more so. Hence, a spatial autoregressive 
model that corrects for spatial error was adopted. Three separate 
spatial error models were run using Geoda 0.9.5-i (Anselin, Syabri, 
& Kho, 2005), with M c „ M cj and M ck as dependent variable, respec¬ 
tively. A queen connectivity structure was used to construct the con¬ 
tiguity weight matrix - that is, zones are considered neighboring if 
their boundaries have at least one point in common. In each of the 
three models, an increase in log likelihood and a decrease for the 
Akaike info criterion was observed, suggesting an improvement of 
fit over the OLS model. The results of the spatial regression models 
are shown in Table 1. 

The spatial autoregression coefficient lambda has a positive ef¬ 
fect and is highly significant in each of the models. When density 
of the working population is controlled for, a higher level of M ci 
is associated with a low job-housing balance and low average com¬ 
muting distance (p < 0.01). The level of Mg, on the other hand, in¬ 
creases with increasing density of the working population, spatial- 
functional diversity and job-housing balance (p<0.01). Surpris¬ 
ingly, proximity to a major city center and distance to a highway 
does not have a significant influence on M a and M g when other fac¬ 
tors of commuting and urban form are taken up in the spatial 


1 The number of departing commuting trips per km 2 is used as a proxy for the 
density of workers in an area. 

2 jhb is calculated as the ratio of the number of arriving commuting trips and the 
number of departing commuting trips. 

3 div of zone i is calculated by means of a Shannon index as div t = P n H P n • 
where N is the number of different land-uses in zone i and p„ is the share of the total 
area of zone i occupied by land-use n. Land-uses were derived from STRUCTNET, a data 
set produced by the Belgian National Geographical Institute (NGI) containing all 
buildings present on the Belgian topographic maps of scale 1:10,000, along with their 
categorization (i.e. industry, school, hospital, public service, etc.). 

4 Major cities in Flanders included the thirteen central cities and Brussels. 
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Table 1 

Results of spatial error regression model. 


Dependent variable: 

M d 



Dependent variable: 

Mcj 
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Mck 



Variable 

Coefficient 

Std. 

error 

v- 

value 

Variable 

Coefficient 

Std. 

error 

P~ 

value 

Variable 

Coefficient 

Std. 

error 

P~ 

value 

Constant 

0.846 

0.064 

0.000 

Constant 

0.629 

0.096 

0.000 

Constant 

1.630 

0.040 

0.000 

dens 

0.775 

0.023 

0.000 

dens 

0.918 

0.036 

0.000 

dens 

0.019 

0.005 

0.000 

div 

-0.055 

0.028 

0.047 

div 

0.173 

0.044 

0.000 

div 

0.010 

0.004 

0.008 

jhb 

-0.128 

0.030 

0.000 

jhb 

0.778 

0.048 

0.000 

jhb 

0.019 

0.004 

0.000 

avgd 

-0.157 

0.037 

0.000 

avgd 

-0.129 

0.058 

0.026 

avgd 

0.003 

0.006 

0.671 

dhigh 

0.071 

0.049 

0.151 

dhigh 

0.087 

0.077 

0.258 

dhigh 

-0.096 

0.009 

0.000 

dcen 

-0.103 

0.062 

0.099 

dcen 

-0.227 

0.094 

0.016 

dcen 

-0.369 

0.024 

0.000 

lambda 

0.575 

0.012 

0.000 

lambda 

0.546 

0.012 

0.000 

lambda 

0.975 

0.001 

0.000 

Pseudo R 2 = 0.554; 




Pseudo R 2 = 0.544; 




Pseudo R 2 = 0.984; 




N = 6668 




N = 6668 




N = 6668 





regression model. This is different in the model with M ck as depen¬ 
dent variable. M ck was found to be positively associated with den¬ 
sity, diversity, job-housing balance, distance to transportation 
infrastructure and proximity of central cities, while the average 
commuting distance did not show a significant effect on M ck 
(p < 0.01). The pseudo R 2 (0.984) for the M ck model is much higher 
than that of M d (0.554) and M CJ (0.544), suggesting a better model 
fit. Also, the relatively high coefficient of dcen indicates that the de¬ 
mand for social interaction is primarily determined by the proxim¬ 
ity of a major city center (see also Fig. 7). 

4.2. Social interaction potential associated with a given commute 

Besides this overall spatial structure of M d , M CJ and M ck , we have 
also explored how M d j may add to our understanding of the geog¬ 
raphy of social interaction potential. Since this metric presumes an 
a priori type of commute from i to j, this exploration can only be 
done anecdotally. Hence, Figs. 8 and 9 show an example of the 
amount of opportunity for social interaction for people living at a 
particular location in East Flanders while working in the inner city 
of Ghent. In other words, the maps represent where workers in the 
inner city should reside to attain a particular level of social interac¬ 
tion potential. It should be noted however that if no commute is 
observed between a zone and the inner city (i.e. Py = 0 in Eq. (8)), 
then the metric attains zero, indicating that there is no contribu¬ 


tion of that zone’s population to the aggregate measurement of 
social interaction potential. For the remaining zones, we have ex¬ 
pressed M C y in two ways: per square kilometer (Fig. 8) and per per¬ 
son (Fig. 9). The former shows a zone’s absolute contribution to M CJ 
of the inner city, while the latter depicts each zone’s contribution 
after controlling for the influence of the size of the working popu¬ 
lation. Generally, social interaction potential decreases with 
increasing distance between the home location and the inner city, 
with the effect of proximity being more articulated after control¬ 
ling for the size of the working population. This means that work¬ 
ers who have a longer commute tend to have less social contact 
options. However, it also appears that even for those having a long¬ 
er commute, living in an urban environment is more beneficial for 
social interaction potential than living in a rural environment. 
Therefore, planning additional residential areas in small urban 
areas may be a valid policy intervention from a social sustainability 
point of view, even if these urban areas have a low job-housing 
balance. 

Zones with a large absolute contribution to the social potential 
of workers in the inner city are primarily encountered in other 
neighborhoods in Ghent surrounding the inner city - the port area 
being the exception - as well as in mainly residential municipali¬ 
ties such as Lovendegem, Evergem, Nazareth and Gavere (see 
Fig. 8); municipalities which are primarily reliant on services and 
jobs in Ghent. This example illustrates that this type of analysis 











T. Neutens et al./Computers, Environment and Urban Systems 41 (2013) 318-331 


327 


# Inner city of Ghent 
-Highway 

'/f- no commute 
^ Waterbody 

Mcij/km 2 

0.0002 - 0.0029 
^ 0.0030 - 0.0081 
0.0082 - 0.0196 

* 0.0190 - 0.0656 

4P 0.0050 - 25.7065 




N 

A 


0 5 10 Km 

I i i i I 



'EekicJ 


ieTdoerlidke] 


F 33rschoot 


[SmtTNikiaas 


Evergem 


■Knes'seiatei 


rZomergem 
kl LOT^idegenr 


Lochristi 


lamm'el 


INevela 


pStellrergen Laam^t^j^B 

l, 

1 e*. Melle^'etter'en’Wiweler 

WSrelbeke W 


BmUWartens-Latem» 

BfiRinti 


iBCiaoenhout 


r KP®ke' 


,Nazareth* 


Gavere 


Zmgen 


t Sint-Lieve ns-Ho utemi 


waaltOTgenderleeuw 


Wortegem-Peiggern' 


M&Wp'riallBBKM 


Ltseraargsirergeci’ 


Fig. 8. Social interaction potential per square km of living in East Flanders for people working in the inner city of Ghent (Mq,/km 2 ). 
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Fig. 9. Social interaction potential per person of living in East Flanders for people working in the inner city of Ghent (M„j/person). 


may for example be useful to inform policy makers’ decision about 
the social sustainability of location alternatives for new residential 
neighborhoods, provided that the work location of prospective 
inhabitants is known. 

Taken together, Figs. 5-9 show that there are marked differ¬ 
ences in the potential to interact with others in Flanders. Depend¬ 
ing on the metric used, the geography of the social potential 


associated with a particular home location, work location or com¬ 
mute pattern can be explored. 

4.3. Social interaction potential within major cities 

Fiaving explored the general spatial structure of the potential 
for social contact, we will now examine the differences of this 
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Fig. 10. Clusters and outliers of social interaction potential of M a - (left), M a — M CJ (middle) and (right) in Brussels, Antwerp, Ghent, Bruges and Leuven. 
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potential within and between major cities in the study area using 
exploratory spatial data analysis (ESDA). In particular, we will use 
Anselin’s (1995) local indicator of spatial autocorrelation (LISA) to 
identify statistically significant clusters and outliers in the geogra¬ 
phy of social interaction potential. The result of this analysis yields 
five cluster types: “high-high”, “low-low”, “low-high”, “high- 
low” and “not significant”. A high-high cluster indicates spatial 
clustering of high values, whereas a low-low cluster indicates spa¬ 
tial clustering of low values. These cluster types reflect positive 
spatial autocorrelation or local spatial clustering of similar values. 
Likewise, a low-high outlier implies that low values of social inter¬ 
action potential are associated with high neighboring values, 
whereas a high-low outlier implies that high values of social inter¬ 
action potential are associated with low neighboring values. These 
outliers are equivalent to negative spatial autocorrelation as they 
point to locations that are different from their neighbors in terms 
of social interaction potential. Finally, the category “not signifi¬ 
cant” refers to the absence of spatial dependence - that is, the spa¬ 
tial pattern of social interaction potential does not differ 
statistically from what could have been expected from a random 
spatial distribution. 

The LISA 5 clusters and outliers of M n , M ci -M c j, and M ck are shown 
in Fig. 10 for Brussels and the four largest cities in Flanders 
(Antwerp, Ghent, Bruges and Leuven). With regard to the social 
potential of a residential location (M a j, low-low clustering is pri¬ 
marily observed in peripheral and port areas. This is for example 
clearly visible for Antwerp and Bruges. High-high clusters of M a 
can primarily be found in the historical centers of the cities (except 
in Brussels), supporting the conjecture that compact living increases 
the social interaction possibilities of a neighborhood. In Brussels 
however, high-high clusters of M c , environ the historical center, 
while the center itself is mainly dominated by low-low clusters of 
M c j-M c j. This is because the historical center contains many neigh¬ 
borhoods where the number of jobs strongly outweighs the number 
of residents (Riguelle, Thomas, & Verhetsel, 2007). This job-housing 
imbalance is, to a smaller degree, also observable in Antwerp and 
Ghent where many low-low clusters of M ci -M c j are encountered in 
the city center. In Ghent, low-low clustering of is addition¬ 

ally present at the intersection of two important highways where 
many employment opportunities are concentrated. In Bruges and 
Leuven, on the other hand, we find low-high outliers at clearly dis¬ 
tinguishable areas in the city. After closer inspection of these areas, it 
appears that these correspond to industrial sites, in the case of Brug¬ 
es, and the university campus and the hospital, in the case of Leuven, 
where the contribution to social potential stems from workers rather 
than inhabitants. Finally, regarding the total amount of social inter¬ 
action opportunity that can take place in the city (M ck ), some well- 
defined pockets of high-high clustering can be observed in the his¬ 
torical centers of the considered cities. These pockets designate areas 
where social interaction intensity is expected to be very high. 

In addition to these internal spatial differences, we can also 
compare cities in terms of their contribution to the overall social 
interaction potential in the study area. In order for the reader to 
get a feel of the spatial structure of the considered cities, Table 2 
summarizes the cities’ contributions to M c , and M q - in the study 
area, the sum of social interaction potential over all TAZs in a city, 
the coefficient of variation 6 in social interaction potential across 
these TAZs, as well as some important city characteristics including 


5 The clusters and outliers have been calculated using an inverse distance 
weighting scheme: all features impact all other features but the farther away, the 
smaller the impact. Other schemes to specify how spatial relationships among 
features are conceptualized were tested, but did result in similar patterns. 

6 The coefficient of variation is a normalized, dimensionless measure of statistical 
dispersion that is calculated as the ratio of the standard deviation to the mean. 
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■ Relative contribution to overall Mcj ■ Relative contribution to overall Mci 

Fig. 11. Contribution of the central cities to the overall social interaction potential 
of the studied region. 



■ Normalized relative contribution to overall Mcj 

■ Normalized relative contribution to overall Mci 


Fig. 12. Contribution of central cities to the overall social interaction potential of 
the studied region relative to the contribution that is expected on the basis of their 
(residential and working) population shares within the region. 

area, population size, road network density, number of origins/desti¬ 
nations of commuting trips, average commuting distance and resi¬ 
dence/employment density gradients. 7 

Fig. 11 shows the contribution to the total M n and M CJ produced 
by each of the central cities. 8 In total, the central cities account for 
21% of M ci and 26% of M CJ - in the study area (Fig. 11). This difference 
tends to suggest that the employment function of the central cities is 
more important than their residential function. However, the social 
interaction potential of workers is more varied than that of residents 
(Table 2) and may in part be explained by the fact that the distribu¬ 
tion of jobs is more centralized than the distribution of residents as 
indicated by the density gradients. The contribution of Brussels to 
the social interaction potential in the region equals over one fifth 
for Md and one third for M cj (Table 2), emphasizing the major role 
of Brussels as employment center, it also implies that the longer 
commute distances to Brussels (Table 2) are strongly inferior to 
the effect of the high residential and working density. Of all central 
cities, Antwerp provides by far the most opportunities for social 
interaction for both workers and residents, followed by Ghent and 
Leuven. Social interaction potential in these cities is supported by 
a dense transport infrastructure. 


7 Following Li (2010), these density gradients have been determined using a 
negative exponential density function: D = Doe where D is the residential/ 
employment density; D 0 is the residential/employment density at the city center; y 
is the density gradient and S is the Euclidean distance to the city center. The larger the 
density gradient, the more centralized the urban spatial distribution of residences/ 
employment is. 

8 For visual purposes, the contribution of Brussels is not shown in Figs. 11 and 12, 
given its outlying value. 


Evidently, a city’s contribution to the overall interaction poten¬ 
tial within the study area depends largely on its population size. 
Therefore, we have divided a city’s relative contribution by the pro¬ 
portion of people living and working in that city for M ci and M CJ , 
respectively. This enables to explore whether the contribution of 
a city to the overall interaction potential in the study area deviates 
from what could have been expected on the basis of its population 
share in the study area. Table 2 shows that the contribution of 
Brussels to M ci is one and a half times as much as what could have 
been expected based on its residential population share. The nor¬ 
malized contribution of to the overall M CJ -, on the other hand, is 
even 2 and a half times larger than expected. Of all central cities, 
Leuven makes the largest contribution to the overall social interac¬ 
tion potential relative to what is expected (Fig. 12), presumably 
owing to its proximity to Brussels. This is also true for Aalst, Ghent, 
Antwerp and Sint-Niklaas. Aalst, which has a very low job-housing 
balance, is the only city for which the normalized contribution to 
M C j exceeds that to M ci . Roeselare, Kortrijk, Hasselt, Genk, Bruges, 
Turnhout and Oostende, for their part, are less sustainable in terms 
of social interaction potential for both living and working. Interest¬ 
ingly, thus, even after controlling for population size, cities within 
or close to the Flemish Diamond still provide more social interac¬ 
tion possibilities than more peripheral cities in the region, not¬ 
withstanding the longer commutes to these cities. 


5. Conclusion 

This paper has shed new light on the contested relationship be¬ 
tween urban spatial structure and social interaction. The study dif¬ 
fers in at least two important aspects from earlier empirical studies 
in this area, which have predominantly concentrated on the anal¬ 
ysis of realized social behavior at the level of residential neighbor¬ 
hoods. First, this study has focused on potential rather than 
realized social behavior. This shift in focus is motivated by the be¬ 
lieve that increasing individuals’ opportunities for social interac¬ 
tion can more easily be achieved through transport and land-use 
policy than the manipulation of their preferences and actual 
choices - a stance that has long since been defended in time- 
geographic work. Second, the method applied in this paper enables 
to study social contact at a much larger scale (i.e. the scale of an 
entire metropolitan region) and accounts for potential social 
behavior outside the residential neighborhood. 

Our findings clearly showed that there is a strong relationship 
between urban spatial structure and the potential for social inter¬ 
action. In general, it was found that the opportunities for after¬ 
work social contact depend on (working) population density, 
diversity, average commuting distance and job-housing balance. 
The demand for social interaction, on the other hand, was found 
to be positively associated with proximity to major city centers, 
distance to transportation infrastructure, density, diversity and 
job-housing balance. Given the centrality of city economics in the 
study area, locations with a significantly high potential for social 
interaction are spatially concentrated in the historical centers of 
major Flemish cities and Brussels. Besides spatial intra-city differ¬ 
ences, we also observed important mutual differences in social 
interaction potential. Cities located near the economic nexus of 
the study area provide more social interaction opportunities for 
both residents and workers, even after controlling for their resi¬ 
dential and working population shares. 

The above findings should however be seen in light of the well- 
acknowledged limitations discussed in Section 3.4. While address¬ 
ing these already constitutes an important first avenue for future 
work, several advancements of the methodology also come to the 
fore. One potential advancement is to refine the metrics proposed 
in this paper by incorporating such aspects as travel time decay, 
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timing preferences (see e.g. Neutens, Schwanen, & Miller, 2010) 
and the spatial distribution of urban facilities for social interaction 
(e.g., restaurants, bars, sports accommodations, etc.). Another is to 
measure social interaction potential at a more disaggregated level. 
This can be realized by deriving individuals’ prism constraints (i.e. 
home/work locations and start/end times at work) from activity- 
travel diaries. By combining this information with aggregate 
commuting flow data, we would be able to calculate the social 
interaction potential of each respondent of the travel diary data 
set, given the realized commuting flows in the area. We can then 
examine whether these individual-based interaction metrics are 
good predictors of the propensity to conduct social activities and 
the locations where these are taking place. Alternatively, one can 
draw on social network data of ego-alter relations (see e.g. 
Carrasco & Miller, 2009). As this type of data generally provides 
the home and work locations of a respondent’s social activity part¬ 
ners, it can be investigated how spatio-temporal constraints affect 
a person’s ability to interact with his/her social network. Finally, 
additional efforts are necessary to compare our findings with those 
in other regional contexts, including non-European cities with 
different morphologies. It is our intention to address these and 
related issues in the near future. 
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